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New expressions are derived for photoemission, which allow experimental electron energy loss data 
to be used for estimating losses in photoemission. The derivation builds on new results for dielectric 
response and mean free paths of strongly correlated systems of two dimensional layers. Numerical 
evaluations are made for BiiSr-zC aCuiO% (Bi2212) by using a parametrized loss function. The 
mean free path for Bi2212 is calculated and found to be substantially larger than obtained by 
Norman et al in a recent paper. The photocurrent is expressed as the convolution of the intrinsic 
approximation for the current from a specific 2D layer with an effective loss function. This effective 
loss function is the same as the photocurrent from a core level stripped of the dipole matrix elements. 
The observed current is the sum of such currents from the first few layers. The correlation within one 
layer is considered as a purely two-dimensional (2D) problem separate from the embedding three 
dimensional (3D) environment. When the contribution to the dielectric response from electrons 
moving in 3D is taken as diagonal in q space, its effect is just to replace bare Coulomb potentials 
in the (3D) coupling between the 2D layers with dynamically screened ones. The photo electron 
from a specific CuO layer is found to excite low energy acoustic plasmon modes due to the coupling 
between the CuO layers. These modes give rise to an asymmetric power law broadening of the 
photo current an isolated two dimensional layer would have given. We define an asymmetry index 
where a contribution from a Luttinger lineshape is additive to the contribution from our broadening 
function. Already the loss effect considered here gives broadening comparable to what is observed 
experimentally. Our theory is not related to the loss mechanism recently discussed by Joynt et al, 
which adds additional broadening beyond what we calculate. A superconductor with a gapped loss 
function is predicted to have a peak-dip-hump lineshape similar to what has been observed, and 
with the same qualitative behavior as predicted in the recent work by Campuzano et al. 
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I. INTRODUCTION 

Photoemission spectroscopy (PES) is an important 
tool to understand the electronic structure of strongly 
correlated quasi two-dimensional systems like high T c su- 
perconductors. Most theoretical work concentrates on 
two dimensional model systems, and when the theoretical 
results are compared with PES the three-dimensionality 
of the actual experimental samples is only schematically, 
if at all, taken into account. Further almost all discus- 
sions are based on the sudden approximation (SA), and 
do not consider extrinsic losses and interference effects. 
For recent work on strongly correlated systems beyond 
SA we refer to Refs. §, §, f§, @. 

We define SA as the bulk one-electron spectral func- 
tion augmented with dipole matrix elements. This ap- 
proximation is exact in the high energy limit for isolated 
systems like atoms and molecules. For solids, where the 
electrons come from a surface region and the mean free 
path is an important feature, SA is never valid, not even 
at high energies. Here the correct high energy limit is 
a convolution of the sudden approximation and the loss 
function (SA*LF). SA is particularly valuable when we 
only look for peak positions such as quasi-particlc ener- 
gies (e g for bandstructure mapping). There are indica- 



tions that also quasi-particle lineshapes are well repre- 
sented. ]5| When it comes to spectral properties over a 
more extended energy region, which is important for e 
g strongly correlated systems, SA can no longer be re- 
lied on. For core level photoemission from weakly cor- 
related systems like metals and valence semiconductors 
SA*LF correctly describes the satellite intensities only in 
the keV region, while the asymmetric quasi-particle line 
shape (in metals) is given correctly by SA already at low 
energies || . For localized strongly correlated systems SA 
is reached rather quickly, say at 5-10 eV above threshold. 



We analyze the three dimensional dielectric response 
of a stack of strongly correlated 2D sheets in the (x, y) 
plane, embedded in a 3D background. We then assume, 
as expressed in Eq. that the response to the total 
electrostatic potential is given by a sum of a 3D part 
and a 2D part. With the 3D part depending only on the 
coordinate difference in 3D, and the 2D part on the differ- 
ence in 2D, the dielectric function is obtained on a closed 
form. This closed form allows us to find an approximate 
relation between the electron energy loss function and 
the dynamically screened potential W . The relation is 
only approximate since energy loss is related to the diag- 
onal part (in q-space) of the dielectric function, while we 
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need the non-diagonal Im W (z, z' . Q, ui) (or equivalently 
Im W (q z , q' z , Q, u>)) for the loss problem in photoemis- 
sion. In PES we need to know lmW(z, z' ,Q,u>) in the 
presence of a surface, while the loss data are obtained 
from a bulk sample. This calls for additional approxima- 
tions. Our numerical evaluations concern Bi2212 and are 
based on a parametrization of the loss function given by 
Norman et al Q. We however include dispersion in the 
dielectric function, which makes the mean free path much 
longer. We use atomic units with |e| = h = m = 1, and 
thus e g energies are in Hartrees (27.2 eV) and lengths 
in Bohr radii (0.529 A). 



£ fc = k 2 /2, u n 



lin (kq-q 2 /2,k 2 /2-k 2 F /2) 



In a solid at lower energies we should use Bloch functions 
and not plane waves for the scattered electron. However 
calculations by Campillo et al show that for copper 
use of plane waves but with a full bandstructure dielectric 
function is a reasonable approximation. In our calcula- 
tions we use Eq. || for the two last terms in the Norman 
et al parametrization. Following Ritchie and Howie ||] 
and many other authors (cf eg ref. we introduce dis- 
persion by replacing uii and Cj in Eq. [I] by Wi(q) and 
Ci(q), 



II. MEAN FREE PATH. 

For the interpretation of photoemission from the 
cuprates the value of the mean free path at energies of 
about 20 eV, where the experiments usually are done, 
is very important. In a recent paper by Norman et al 
j| very short values of the order 2-3 A were obtained 
for Bi2212. Norman et al however neglected the q de- 
pendence in the loss function. In the electron gas case 
neglect of dispersion makes the mean free path about half 
the value with dispersion. When we introduce dispersion 
for Bi2212 we find an even larger effect on the mean free 
path. 

Norman et al used a parametrization of the energy loss 
data on Bi 2 Sr 2 C 'aCu 2 Os (Bi2212) obtained by Niicker 
et al 01, 



uii(q) =0Ji + —, Ci(q) 



2,3. 



We have put kp = for simplicity, which gives a slight 
underestimate of the mean free path. 

The expression for the mean free path in a layered ma- 
terial is 
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where k z and K etc are components perpendicular and 
parallel to the layers, k = (k z , K), k = |k|, K = |K| etc, 
and 
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with parameters (energies in eV) given below 



i d w, Ti 

1 0.164 1.1 0.7 

2 0.476 18.5 13.6 

3 0.345 32.8 17.0 

The first peak at about 1 eV is associated with 2D plas- 
mon excitations, while the large double peak comes from 
essentially 3D excitations since it is similar to what is ob- 
served in Cu metal (c f Ref. Q). The linear rise for small 
lo comes from acoustic plasmons (due to the coupling of 
the 2D plasmons in the different layers), and also to some 
extent from electron-hole excitations. Phonons and other 
low energy excitations cannot be seen in Niicker et al's 
H data since the broadening is too large (150 meV). 

For an electron gas we have the well-known relation 
between the mean free path (sk) and the inverse di- 
electric function Ime -1 (q,ui), 



^3D (sk) 
with 



2 
Trk 2 



dq 
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£ (q,ui) 



du, (2) 



For simplicity we have taken free electron energies. Fur- 
ther considering propagation perpendicular to the layers 
we have K = 0, and no dependence on the angle <p be- 
tween K and Q, 
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In ref. |l0| there is a detailed discussion of the 1 eV fea- 
ture (the first peak). It has a Q 2 dispersion which is a 
signature of coupled particle-holes (plasmons) . They also 
estimate the coefficient theoretically with reasonable re- 
sults. From Niicker et al Q it is clear that the first peak, 
besides dispersing as Q 2 , quickly broadens when Q in- 
creases. For the first term in Eq. we use 



co 1 (Q)=Lo 1 +aQ 2 , r 1 (Q)=T 1 1 



ci (0) 



<4(Q) 



, Qq = 0.13 ait, a = 0.6. 
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We have put [i = 0. With a finite fi, and thus a fi- 
nite fewe should replace y/k 2 — Q 2 with \Jk' F — Q 2 for 
Q < kp in the limits above. Such a replacement makes 
I/A2D smaller and our approximation thus again slightly 
underestimates the mean free path. 
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FIG. 1. The contributions to the inverse mean free path 
1/A from 2D (full drawn) and 3D (dashed) terms in the case 
of Bi2212. 
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FIG. 2. The full drawn curve gives the mean free path A 
from our calculations which include dispersion in the dielectric 
constant, and the dashed curve the results without dispersion 
given by Norman et al [4]. 

In Fig. 1 we show results for the inverse mean free 
paths for the 2D and 3D contributions. It is remarkable 
that the 2D effects, while peaked at 1.5 eV, extend quite 
far, out to some 30 eV. The 3D contribution starts dom- 
inating at about 10 eV. The mean free path A, given by 
1/A = 1/ A 2 £> + 1 /A30 , is shown in Fig. 2. The maximum 
in I/X3D is reached at about 100 eV, where the mean free 
path has its minimum of some 5 A. It seems to be a uni- 
versal feature that the minimum mean free path is about 
5 A at an energy about 3-4 times the energy where the 
loss processes are strongest, as can be seen from tabu- 



lations of loss functions [gdj and mean free paths. |12| 
The qualitative behavior of the 2D and 3D contributions 
in Fig. 1 are similar to what has been obtained in RPA 
calculations for the layered electron gas. JT]| We remark 
that Norman et al besides the inverse mean free path also 
calculated the background in PES in the traditional way 
from the extrinsic losses only, following the common con- 
vention to take the background as zero at the bottom of 
the main band. If we do a similar background calculation 
the results are very close to the Norman et al's since in 
such a calculation only the shape and not the strength 
of the loss function enters. We emphasize that in our 
treatment of PES later in this paper we include besides 
the extrinsic losses also the intrinsic ones and the inter- 
ference terms, which gives a radically different result for 
small energy losses. 

III. PHOTOEMISSION 

We are interested in photoemission from the CuO lay- 
ers. The layers are regarded as localized systems embed- 
ded in a 3D environment. The crystal surface is taken to 
be parallel with the layers and in the (x, y) plane. We 
take the electrons in the 2D layers as separate from the 
other electrons, and write the state vectors for the initial 
and final states as 

\N B ) \N 2D ) , \N%, ax) \N 2D - 1, s 2 ) |k) . (4) 

Here \N 2 d) is the state vector for the electrons in one 
particular layer at a distance zo (zq > 0) from the sur- 
face (the one from which the photoelectron comes), and 
\Nb) the state vector for the remaining (bulk) electrons 
which move in 3D. \N 2 d — 1, s 2 ) is an excited state s 2 of 
the particular layer, and \Ng,si) an excited state si of 
the bulk electrons. The star indicates that these electrons 
move in the presence of a localized hole at r = (0, 20)1 
and thus is an cigenfunction of a different Hamiltonian 
than that for |Ab). Finally |k) is the photo electron 
state. One may argue that the hole should be extended 
over the 2D layer rather than sit in one point. However 
even in weakly correlated solids correlation effects give 
rise to satellites corresponding rather to the removal of 
the electron from a point than from an extended region. 
Q For the strongly correlated systems considered here 
the band is quite narrow and thus the atomic functions 
building the Bloch functions have small overlaps, which 
makes our approximation of a localized hole even better. 
We consider processes when the photo electron energy is 
high enough that we have reached the sudden limit as far 
as the excited layer is concerned (about 10 eV according 
to Ref. §). 

The expression for the PES transition amplitude then 
becomes |D| 

r(k, Sl ,s 2 ) =^2(k\(N B , Sl \(N 2D -l,s 2 \ (5) 
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1 + V 



E- H 



Cl \N 2D )\N B )A\t) 



where \i) and |k) are one-electron states. The state 
\i) = \K.i) \cj>o) is the product of a 2D Bloch state |K;) 
with momentum Kj , and the bound state </>o (z) for 
the motion in the z direction, which only will enter as 
\(j)Q (z)\ = w(z). The operator Cj destroys an electron 
in state "i" , and c; |iV2D} is regarded as a localized state 
concerning its influence on the 3D states. The optical 
transition operator is A, and V is the potential for the 
interaction between the photo electron and the solid, the 
potential that causes external losses. Since V contains 
both operators acting on the photo electron and on the 
solid the expectation value (N B ,s\ \..V..\Nb) is a one- 
electron operator acting on photo electron states. The 
state A \i) generated by optical excitation is considered 
a photo electron state. H is the full Hamiltonian includ- 
ing V, and E is the total energy 



E = E (N 2D ) + E (N B ) + 

= E(N B ,s 1 )+E(N 2D -l,s 2 )+e lc , 



(6) 



where ujphot is the photon energy and e k = k 2 /2 the 
photo electron energy. The photocurrent is proportional 
to 

Jk (oJphot) = ^ \T (k, Si, S 2 )| 2 S (ujphot + E S2 - UJ Sl - £ k ) . 



with 



E(N B , Sl )-E(N B ). 



E{N 2D )-E(N 2D -l,s 2 ) 



For the Hamiltonian H we take, 

H = H 2D + H QB + h + V, 

where H 2B . describes the pertinent 2D layer, Hq B the 
3D electrons in a quasi-boson representation, h (a one- 
electron operator) the photo electron, and V the inter- 
action between the photo electron and the 3D system 
(the interaction with the 2D system is neglected since we 
assume the sudden limit to apply here). Explicitly we 
have, 



H, 



QB 



w.aLa. 



V h P h , V h = Y,Vh 



4) 



potential from the hole in the 2D system. The functions 
V s (r) are fluctuation potentials, discussed at length in 
Refs. d and H. 

Say that we somehow can calculate the photo current 
( z 0i ujphot) from one isolated two-dimensional layer 
at a distance zq from the surface, and want to estimate 
the current from this layer when a set of such layers to- 
gether with additional electrons of 3D character form a 
three dimensional crystal. We have to account for the 
shake up in the 3D surrounding of the layer as well as 
the losses the photo electron can have on its way out to 
the surface. In Appendix A we show that the photocur- 
rent then can be written as a convolution between the 2D 
current (zqAjj) and an effective broadening function 
P k Oo, w) (Eqs ||and |H), 



2D 
k 



(z , lu') F k (2 , LOphot ~ w') dJ. 



(7) 



A delta function peak 5 (u> — Eq — e k ) in J^ D (zq, uj) will 
hence give a contribution P k (zo,u>phot — £0 — e k ) to the 
photocurrent. In core electron photoemission we have a 
similar expression with J^ D (zq,uj) replaced by the ex- 
pression for the current from a core level in an isolated 
ion (essentially a delta- function) . 

The common EDC experiment gives the current 
for fixed photon energy u> P hot as a function of the 
electron energy e k for a given direction of k (or 
K, where k=fc z z+K). We are thus interested in 
P k (20, Wphot — £0 — £k) as function of e k in the range 
Wphot — £0 > £k > 0. Since Pk (zo,oj) varies fairly slowly 
with k for fixed uj, Pk (20, uj) as function of uj for fixed k 
describes the photo emission (in a limited energy range) . 
We will mainly discuss the properties of P k (2 ,w) as 
function of uj for fixed k. 

The effective broadening function to second order in 
V s is found to be (Appendix A), 



P k (z ,oj) = e 



6(uj) 



a(k, z ;uj) 



uj 



(8) 



where 
a(k, z ;uj) 



uJ 



E 



/ (k, Q, uj, 2 ; 2) V (q z ,Q>z) dz 



xS (uj — uj s ) . 



(9) 



V = E ^kk4ck' (a. + a\) , - (k \V S \ k'> , 



skk' 



/ (k, Q, uj, z ; z) 



8(z- zq) '_ ik.z 

uj in 



(10) 



V,t =V s (0,z ). 

Ph is a projection operator which gives 1 for states with 
a hole in the 2D system, and otherwise, and Vh is the 



This expression is the same as in Eqs. 26, 27 for the core 
electron current in Ref. ||. The function V (q z ,Q,z) in 
Eq. [5] is the fluctuation potential giving the coupling 
between the photo electron and a density fluctuation 
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s = (q z , Q) with energy u a . In / (z) the first term gives 
the intrinsic or shake up contribution to the amplitude, 
while the second term gives the contribution from losses 
when the electron propagates from the layer at zq > 
to the surface at z — 0. The quantities k and k z are the 
(complex) momenta in the z direction of the photo elec- 
tron when inside the solid before and after it excited the 
density fluctuation s having parallel momentum Q. The 
photo electron momentum outside the solid is fc 2 z+K 
and its energy £k = {k 2 z + K 2 ) /2. Further Vq is the (neg- 
ative) inner potential, and Ti and T2 arc the dampings 
before and after emitting the excitation s. It is easy to 
derive expressions where the plane waves e lkz and e~ lKZ 
are replaced by (damped) Bloch functions, and doable to 
find expressions where the bandstructure also is present 
in the lateral motion. 

Since (u>) is quadratic in the fluctuation potentials 
V s , we can relate it to the dynamically screened poten- 
tial W . For the imaginary part of W we have (c f Eq. 49 
in Ref. |), 

ImW(z,z';H,R,';uj) = -ttJ^ V s (r)V s (r')S(uj - u a ). 

S 

With ^ s (r) = e iQR V (q z , Q,z) 

Im W(z, z'; Q; u>) = -irA^V (q z , Q,z) V (q z , Q,z') 

xS(uj - uj s ), 

where A is the normalization area of the planes. In an 
exact treatment the V (g z ,Q,z) can be chosen real, and 
we see that Im W(z, z'; Q; lo) is symmetric in z and z' . 
Comparison with Eqs. || and ^ shows that 



a(k, z ;u) 



—y 

IT A ^ 



Q 



/(k, Q,w, z ;z) 



x ImPy"(z, z'; Q; uj)f (k, Q, w, z ; z')* dzdz'. 

To simplify the calculations we relate Im W to the mea- 
sured loss function. The loss function however is con- 
nected with losses in the bulk, and we also have to find an 
approximate relation between lmW bulk and ImW 8 '"''. 
This was done in Ref. H by using the Inglesfield simpli- 
fied expression for the fluctuation potential, 



F{q z ,Q,z)F(q z ,Q,z') 



2tt _ _ 

xlmW^ lk {q z ,Q,ui)dq z , 



n 



where 

F (q z , Q,z) = 2 [cos (q z z + q ) 



cos q e Qz ] 



6 q = arctan — . 



This means that for the strength of the coupling we keep 
the bulk expression, while for the spacial part we have a 
bulk function (here plane wave) which is modified to be 
zero at the surface. The relation to the loss function is 



ImW^f {q z ,Q,")=v{q z ,Q)lm 



1 



e? D (g z ;Q;^)' 



For 2D excitations we can do a similar modification of 
the bulk fluctuation potential to make it zero at the sur- 
face. When we take w (q z ) = 1 we have (see Appendix 
B) 

ImW™ rf (z,z>;Q;u) 

= T- I' C yT (««' - z o) yT ^ - z o) 
27r Jo 

T xo(^;Q;w) 

x Im dq z 

c 

for the contribution from the layer at zo, where 
Imx(q; Q; w) is related to the loss function by 



Im 



Xo(g z ;Q;w) 



1 



Im ■ 



-1 



2v(q z ,Q) £ 2D (gz . Q;w) > 
and V (q z , Q,z) is the fluctuation potential, 
V r (q z ,Q,z) = 2Be[V (q z , Q; z) exp {-iq z z + i<f> (z ))} . 

(11) 

Here V p (q z ,Q;z) is a well-known periodic potential 
(V'P = exp [iq z z) V, with V defined in Eq. ||) 



V p (q z ,Q;z) =J2v{q z + G,Q)w(q z + G), 



-iGz 



(12) 



_ 2irce lq * z sinhQ (c - z) + e -^ c smh Qz 
Q cosh Qc ~~ cos q z c 

The explicit expression follows when the form factor 
w (Qz) ( c f Appendix B) is taken as 1, and is valid only 
for < z < c. The phase <fi (z ) in Eq. [TT] is chosen to 
make V r zero at the surface, V r (q z , Q,zq) = 0. 
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FIG. 3. The periodic functions Re V p (z) and Im V p (z) for 
z/c in the interval (0, 1), where c is the lattice constant and 
V p is defined in Eq. Il2l We have taken some typical values, 
q = 0.03 and Q = 0.1 (for Bi2212 c = 29.1 and a = 10.2 
which gives 7r/c = 0.11 and ir/a — 0.31). 

In Fig. 3 we plot Re V p (z) and Im V p (z) for some typ- 
ical values of q z and Q, and in Fig 4 we show V r (z) for 
the same parameter values. The sharp peak with a singu- 
lar derivative in Re V p (z) at z = is smoothed if we take 
uu (q) 1 . For a typical binding energy of 3 eV and an ex- 
ponential wave function, we have w (q z ) — a 2 / (a 2 + q 2 j, 
with a = 0.9. Typical values of q z and Q are n/c and 
it /a. The lattice parameters for Bi2212 are c = 29.1 
and a = 10.22 which makes ir/c — 0.11 and it /a = 0.31. 
Jr6| We can also compare with the cut-off parameter for 
the collective excitations in Bi2212 discussed in the sec- 
tion on mean free path, Qq = 0.13. Thus a is substan- 
tially larger than q and Q, and it is hence reasonable to 
take w (q z ) = 1. We note that the values of KeV r (z) 
at the first two Cu layers are substantially smaller than 
the maximum value of He2V p (z). An approximation 
with bulk potentials cut at the surface clearly can give 
very large and spurious effects unless we go to so ex- 
tremely high energies that the mean free path becomes 
much larger than the lattice parameter c. 
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FIG. 4. The fluctuation potential V r (q, Q;z- z t ) in Eq. [y] 
for z/c in the interval (0,2.5) and for Zi = 0.39c and 0.61c, 
the distances from the surface of the first two CuO layers. 
The potentials are zero at the surface and the cusps come at 
CuO layers. The maximum possible value of \V r \ is 2V P (0). 

When there are no low energy excitations, like for an 
insulator or for a metal when the electron-hole excitations 
are neglected and only the plasmons are kept, the over- 
lap between the initial ground state and the completely 
relaxed ground state in the presence of a localized hole 
potential is finite. In a quasi-boson treatment we have 



\(n* b mnbM 



V s (z ) 



A partial summation of the perturbation expansion in V s 
(or a cumulant expansion) gives 0, 



dt 
2^ 



e- iut exp 



a (k, z ;uj') ■ 



-du' 



(13) 



This expression correctly reproduces the edge singularity, 
l/ti> 1_Q ^ k ' z ° ;0 ), and also gives the second order satellite 
term in Eq. |^. In the high energy limit and the plasmon 
pole approximation (the electron-hole part is then not 
included) it can be shown analytically that || , 



duj = a + zq / X, 



(14) 



and Eq. [13] thus also gives the correct prefactor in this 
limit. 
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FIG. 5. Results for «2D (w) and (xzd (w) (c f Eqs. and 
|lg| ). The parameter values are obtained from energy loss 
data for Bi2212. The curves in figs. 5a,b give contributions 
to ct2D (w). The symbol mn refers to a contribution when 
the fluctuation potential is centered at layer m (and m+2, 
m+4 etc) and the photocurrent comes from layer n. Thus the 
dashed curve in fig. 5a (m=l, n=2) refers to a contribution 
from the fluctuation potential centered on the layer closest 
to the surface when the photocurrent comes from the second 
layer. Fig. 5c shows the total contributions to 0215 (w) when 
the current comes from layers 1 to 4, and fig. 5d the corre- 
sponding contributions to a^D (w). The curves in the inset 
are ad hoc adjusted to take out the unphysical low energy 
part coming from a schematic parametrization. The photon 
energy is 1 au. 

With an electron-hole continuum, Eq. |l4] no longer 
can hold since the integral diverges. We then split 
a (k, u>phot) k>) m a 2D part from the excitations in the 



G 



layers, and a 3D part from the remaining excitations, 
a = ct2D+&3D- The 3D contributions in Eq. [l] have been 
smoothly deformed to be zero for to < u t h =0.1 since the 
metallic excitations come from the layers. To give a good 
representation of the experimental loss function this de- 
formation should be compensated by a small increase in 
the 2D term, but this is a minor effect which we have 
omitted. Now the integral J u dioa^o (k, u> p h t;uj) /to 
converges, and we have checked numerically that in the 
high energy limit 



a 3D (k, z ;u) 



doj = a 3 % r 



(15) 



Eq. |l6| guarantees the correct dependence on the dis- 
tance z to the layer. While exp (—a 3 n p r ) may not give 
an accurate scale factor, this is of minor importance since 
it does not affect the ratio between the threshold peak 
and the satellite structure. The a 3 n ^ r values depend only 
weakly on zq, and for the first four layers the values are 
0.243, 0.252, 0.260, and 0.261. Collecting our results we 
have 



a 2 D (k, z ; lu) 



7f/c /■ 

3 / <^z / d Q 

tt(2tt) Jo J 
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(17) 



lmx(g 2 ,Q,w), 



where a 3 n p r contains only the intrinsic part, 
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a 3D (k, z ; uj) 
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ImWb (q z ,Q,v) 



The approach to the high energy limit is quite slow (of 
the order of keV), and in our estimates for Bi2212 we 
adopt the expression 





intr 
a 3D 




poo 


x exp 






Jo 






+ 


J ' 



dt 
— < 
2tt 



(16) 



a 2 D (k, z ; u') : dtu' 

lu' 

+ 1 a 3D (k, z q ;l>') — — dw' 
1 lu' 



z0=0.39c - 
z0=0.61c 
z0=1 .39c 
z0=1.61c 



E(eV) 



a 

C\J 

8 t 













full - 




intrinsic 




l\<0.3 








z0=1.61c 




\ \t0.2 


z0=1 .39c 




\toi 


z0=0.61c 




z0=0.39c 



a 

CM 0.3 

0.2 



E(eV) 



3.61c 
3.39c 
2.61c 
2.39c 
1.61c 
1.39c 
0.61c 
0.39c 



E(eV) 




10 20 30 40 50 60 70 

E(eV) 



FIG. 6. The left part shows a comparison of a2D (lu) for 
the pure intrinsic case with the full expression including the 
extrinsic amplitude. The right part shows how the pure in- 
trinsic contributions to c*2d (lu) and a$D (lu) converge towards 
their bulk values. The photon energy is 1 au. 




E(eV) 

FIG. 7. The effective loss function (lu) including both 
2D and 3D contributions (c f Eq. [TJ|). The contributions 
from the different layers are displayed separately. The curves 
are convoluted with Lorentzians, in Fig. 7a with r = 10 meV, 
and in Fig. 7b with T = 300 meV. 

In Fig. 5 we plot results for different contributions to 
the a functions in Eq. [l(| The general shape of the 
CX-2D (w) functions is similar to the electron gas case with 
a flat portion for small lu followed by a plasmon peak 
(c f Ref. flTi] , pp 663-667) . However the magnitudes are 
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different, aiD (0) is fairly large (0.25 — 0.30) compared 
to metals while the plasmon peak is much smaller and 
broader. When we change the parametrization to make 
the 3D terms start at 0.1 au the a.2D (0) values will in- 
crease by some 10%. The olid functions have only a weak 
dependence on photon energy, while the a 3 o curves have 
a much larger dependence. All curves in Fig. 5 are for 
the same photon energy of 1 au. 



> 
H 

CO 

z 

LU 



15 



10 



>- 

H 

CO 

z 

LU 



0.6 



0.4 



0.2 





(a) 




full solution 


\ 


rectangular model 

^T^^^h= 1 | | 


1 2 


3 4 5 6 7 6 




E(eV) 




(b) - 




full solution 




rectangular model 



E(eV) 

FIG. 8. The effective loss function P k (uj) obtained with 
only the 2D contribution from the second layer. The full 
drawn curve in fig 8a shows the full solution obtained with 
the "2" curve in Fig. 5c, while the dashed curve shows the 
result using the rectangular approximation for Q2D (lo) with 
ujo = 0.08 and ao = 0.255 (see text at Eq. |l^). In fig. 8b we 
show o; (1 ~ Q(0)) P k (w). The photon energy is 1 au. 

The left part of Fig. 6 shows the dominance of the 
intrinsic contributions to U2D- As expected the contri- 
butions from the extrinsic terms are larger for the layers 
further away from the surface. The right part shows the 
approach towards the bulk value of the intrinsic contri- 
butions for a.2D and a^o- This approach is considerably 
slower in the 2D case as might be expected from the be- 
havior of the fluctuation potentials (c f Fig. 4) . Compar- 
ing the intrinsic a^D in Fig. 6 with the full a^D in Fig. 5, 
we see that in the 3D case the extrinsic effects dominate. 



full 
3D 



a 3D T i s roughly proportional to 



The difference a 
zq which follows the trend in the high energy sum rule, 
Eq. [L5l In Fig. 7 the contributions from Pk (zq,lo) in 
Eq. Hfl from the first four CuO layers are shown. It is 



clear that most of the asymmetry comes from the layers 
in the first unit cell. The alpha function for the first cop- 
per layer is quite small (Fig. 6), but when the mean free 
path effects arc taken into account, Fig. 7 shows that 
the broadening contributions from the first and second 
layers are comparable. Fig 7b shows an extended energy 
region to illustrate the relative importance of the 2D and 
3D contributions. The integral effect of the 3D contribu- 
tions is much larger, but the peaks in the loss function 
are smoothed out and the 3D contribution is featureless. 
At higher energies we of course also have contributions 
to the photo current from other states than the quasi 2D 
ones in the Copper layers discussed in this paper. 

We now give a qualitative discussion of the effective 
broadening function Pk (<■<-> ) in Eq. [l6| Since the 2D and 
3D contributions add in an exponent we can write Pk (lo) 
as a convolution, 



Pk (z ,to) = e 



-z /\-a' 3 ^ r 



~Pl° (lo')cIlo'. 



For P^ D we make a Taylor expansion, and keep only the 
first term, P^ D (ui) = S (ui) + a 3 o (w) /lo. We have then 
omitted the multiple quasi-boson excitations starting at 
lo = lLOth- Since P£ D is normalized to unity, and consists 
of a peak that is sharp compared to a^o , we can write 



Pk {zq,lo) ~ e 



-z /x— <45) r 



Pl° (z Q ,Lo) + 



a 3D (k, z ;lo) 



To numerically evaluate Pk D (w) we used the integral 
equation loP£ d (lo) — J" cL/a 2 D {lo') P£ d (lo — lo') which 
is easier then to evaluate the exponential expression in 
Eq. |l6|. If we approximate aiD (lo) by a rectangu- 
lar function, o^d ( a - 1 ) = a o@ (<^o — an d broaden with 
a Lorcntzian of width T (FWHM = 2T), we have for 
lo < loq the Doniach-Sunjic expression p8[, 



P k 215 M = C(ao) 



C(a ) 



cos [7rao/2 — (1 — ao) arctan (lo/T)} 
(l + (lo/T) 2 



-■ya 



(a a - ly.LO^T 1 -^ sin [ira ] 



(19) 



where 7 = 0.577 is the Euler constant. The coeffi- 
cient C (a ) in Eq. |l^ was derived in Ref. (see 
Eq. 162). For lo > lo q P£ d (lo) only has a weak tail 
with less than 10% of the norm (for ao < 0.4). Let 
be the lo value for which Pk D (lo) has its maximum, 
and lo\ and L02 the values where it takes half its maxi- 
mum value. We define an asymmetry index j as (ao) = 
(L02 ) / (wmax — wi). An approximate expression 

for j as is 



las (a ) = 



L0 2 



L0 n 



■ LOl 



= 1 + 0.79a 



14.54ag. (20) 
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p2D 



function that broadens a (5-function 

2D 



(lo) is the 

peak in J^ D (ui). If j£ v (lo) has a Doniach-Sunjic sin- 
gular shape the broadening with P£ D (lo) still gives Eq. 
|l9| but with an ao that is the sum of the alphas in J^ D 
and in P£ D (lo). This is so because the time transform of 
a power law singularity t a ° , and a convolu- 

tion in frequency space is a product in time space. 
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FIG. 9. The effective loss function Pk (u>) convoluted with 
Lorentzians of different widths T (5, 10 and 20 meV). All the 



2D contributions P k (ui) in 



Eq. ^ 



from the first four lay- 



a$g r ). The 3D 



ers are summed weighted with exp(— zo/\ 
terms are not included. In the first three panels we have used 
the rectangular approximation, while in the last panel the full 
evaluation from Eq. [TJ] was done. Also the Lorentzians are 
shown. The photon energy is 1 au. 

The a2D functions in Fig. 5 show clear peaks due 
to the plasmon excitations. The peaks are however not 
strong enough to give more than a small bump in the 
P£ D functions. This is illustrated in Fig. 8a which 
shows the full P? D curve and the rectangular approxi- 
mation in Eq. |l9| using ao = 0.255 and luo — 0.08 au — 
2177 meV. In Fig 8b the rectangular approximation 
is illustrated by taking out the singularity and plotting 
P£ D (lo) uj^^oio)) _ The simple rectangular approxima- 
tion without plasmon peak should be useful as a guide 
when other broadening effects are at work. In Fig. 9 we 
show the sum for the first four layers of the 2D contribu- 
tions exp (— z /A — a^ r ) P£ D (lo) broadened with differ- 
ent Lorentzians. The 3D terms are not included except 



for the (all important) mean free path factor. In the first 
three panels with a limited energy region (up to 500 meV) 
we have used the rectangular approximation for the dif- 
ferent «2D contributions. In the last panel with a larger 
energy range the full evaluation from Eq. ^| was done 
since it is superior to the rectangular model for energies 
above 0.5 eV (see Fig. 8). The numerical accuracy at the 
peak is however lower in the full calculation. Also the 
Lorentzians are shown to ease the estimate of the size of 
the asymmetries. It is clear that we have a sizeable line 
asymmetry, and also a long tail extending over several 
eV. The artificial step in the rectangular approximation 
at about 3 eV (Fig. 8) is of little consequence since the 
intensity is small at this energy. The asymmetry index is 
slightly dependent on the Lorentzian broadening T since 
we have summed contributions from different layers with 
different a values. The index is about 2.6 which accord- 
ing to Eq. ^ corresponds to an effective a of about 0.3. 

In the superconducting state the loss function should 
have a gap. We mimic this gap by using a rectangular 
alpha function 



a2D M = a 9(uj - lo sc )6 (luq - lo) 



(21) 



using the same values for ao and loq as in Fig. 8. For 
the gap lo sc we take lo sc = 70 meV. In Fig. 10 we show 
the corresponding P£ D (lo) broadened with a Lorentzian 
of width r = 15 meV. Our choice of parameters is 
only made to illustrate the qualitative behavior to be 
expected. The curve clearly shows the peak-dip-hump 
lincshape found experimentally (for a recent reference see 

eg H) 
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FIG. 10. The loss function P£ D for a gapped spectrum us- 
ing the simple parametrization in Eq. Ell The Lorentzian 
broadening is V = 15meV. 

Recently it has been possible to obtain very accurate 
tunneling data from Bi2212, and it is of interest to com- 
pare these data with the PES satellites (ref. p0|), since 
the tunneling data also show peak-dip- hump structures. 
pl[ PES and tunneling are basically different spectro- 
scopies. There can however be qualitative similarities 







since in both cases the electrons couple to 3D quasi- 
boson excitations like phonons, electron-hole pairs, plas- 
mons, magnons etc. In our treatment of PES we take 
the states of a particular 2D layer as given and study the 
effect to low order of the sudden appearance of a hole in 
the 2D system on the quasi-bosons (intrinsic excitations) 
as well as of the coupling of the photoelectron leaving 
this layer to the quasi-bosons (extrinsic excitations), and 
their interference (c f Eqs. 8-10, or equivalently Eq. 26 in 
Appendix A). We found that the intrinsic contributions 
dominate for small excitation energies. 

Tunneling is traditionally described by the spectral 
function which involves matrix elements of the electron 
annihilation operator between the initial state and the 
excited states. | f22[ , p3| The excited states consist of a 
2D layer state with a hole, and some state of the quasi- 
bosons in the presence of a localized hole. In lowest order 
perturbation theory the probability for a final state with 
excited quasi-bosons is given by the first term in Eq. 10. 
This means that the intrinsic contribution to PES and 
the tunneling currents are the same except for the mean 
free path effect in PES shown in Eq. 8, and the sum- 
mation over momenta in tunneling (giving the density 
of states, DOS) contra momentum conservation from the 
dipole matrix element in PES. As mentioned above, we 
modify this analysis valid for the normal state, by simply 
assuming that the loss function should have a gap in the 
superconducting state. 

In Bi2212 we have a van Hove singularity (VHS) at the 
Fermi level, which makes the difference between DOS and 
momentum conservation of less importance (there might 
actually even be two VHS if the two CuO planes at 3 A 
apart produce a significant splitting). More important 
is that in PES the electrons come from a thin surface 
region (of the order of the mean free path) while in tun- 
neling they may come from an extended region which can 
be hundreds of A, and that the coupling functions V(z) 
are strongly z-dependent. Additionally there are two en- 
ergy gaps (superconducting gap and pseudogap), which 
further complicates the picture. There is thus no way 
that PES and tunneling structures can be quantitatively 
the same, but since the same quasi-bosons are involved, 
there may well be qualitative similarities even though the 
coupling strengths can be quite different. It should also 
be noted that we take the spectral function for the 2D 
system as a sharp peak (the function often calculated by 
theoreticians using say a t— J model), and have no means 
to estimate the relative strengths of the true 2D spectral 
function relative to the loss structure analyzed here. 

In our analysis we have only treated the plasmons for 
the simple reason that the experimental loss data at hand 
did not have resolution enough to show phonons and 
other low energy excitations. If such (q = 0) data ap- 
pear showing additional quasi-bosons one is faced with 
introducing reasonable dispersions, and finding reason- 
able extrapolations of the bulk coupling function to ac- 
count for the presence of the surface. 



IV. CONCLUDING REMARKS 



This paper is concerned with effects of external losses 
in photo emission, and the extent to which the com- 
monly used sudden approximation works for strongly cor- 
related layered materials. We have earlier found that for 
a strongly correlated localized system the sudden approx- 
imation is reached rather quickly, at about 10 eV [GJ. For 
a weakly correlated system on the other hand, like an sp- 
metal or semiconductor, the sudden limit is approached 
very slowly, on the keV scale. [|| The slow approach is 
connected with strong destructive interference between 
the intrinsic and extrinsic mechanisms for plasmon pro- 
duction. The cancellation is particularly strong for small 
momentum plasmons where the long-wave plasmons are 
excited by the average potential from the core hole and 
photo electron, which is zero [ gjj . The asymmetric line- 
shape in core electron photo emission from metals is, on 
the other hand, hardly affected by the external loss pro- 
cesses ||. 

We are interested in energies where the sudden limit 
is reached for the strongly correlated layer from which 
the photo electron comes, and derive an expression for 
the photo current as a convolution of the sudden approx- 
imation for the current from the layer with an effective 
loss function, Pk(uj) (Eq. We assume, as far as the 
loss properties are concerned, that the photo electron 
comes from a localized position. In our specific example, 
Bi22l2, the c value is 15.4 A (neglecting crystallographic 
shear), and almost all contributions come from the first 
unit cell. The two first CuO layers are at 0.39 c and 
0.61 c from the surface (which is between two BiO lay- 
ers). 1 25 With a photon energy of 1 au the maximum 
electron energy (inside the solid) is 1.15 au = 31 eV if 
we take the bandwidth as 0.15 au. Our energy loss cal- 
culations give A = 17.8 a (Fig. 2). The exp(— z /A) 
factor then is 0.53, 0.37, 0.10, and 0.07 for the first four 
CuO layers. We thus expect large photo emission contri- 
butions only from the first unit cell. 

To obtain Pk (w) we use a previously developed method 
based on a quasi boson model, where the electron-boson 
coupling is given by fluctuation potentials related to the 
dielectric response function Jl5|] . We find that the energy 
loss function, which we take from experimental data, can 
be related to the screened potential which we need to cal- 
culate the (intrinsic and extrinsic) losses in photo emis- 
sion. The fluctuation potentials related to the electrons 
in the layers are universal functions, which are easily cal- 
culated (Eqs [y] and [l2]) . They have some resemblance 
to a surface plasmon potential, but penetrate the whole 
solid and have the Bloch wave symmetry. We use the 
real part (or equivalently the imaginary part) of a phase 
shifted bulk potential to get a potential which is zero at 
the surface, and mimics the potential we have in a finite 
solid. The fluctuation potential is integrated over z to- 
gether with a propagation function / (z) (Eq. O) that 
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takes the photo electron out of the solid. This integral 
is in turn integrated with the loss function (taken from 
experiment), Ime -1 (q, Q, cj), to give functions o^d (w) 
and a^D (w) which are simply related to the effective loss 
function Pk (<^) (Eq. |l6|). When we use plane waves 
instead of Bloch functions in the propagation function 
f{z), all specific materials properties are embodied in 
the loss function. The propagation function has both an 
intrinsic and an extrinsic contribution which interfere. 



From Eq. 16 we see that Pk (^) is scaled down with zq, 
the distance of the layer from the surface, while the fluc- 
tuation potentials increase with z. The reason for that 
increase is that the boundary condition forces the fluc- 
tuation potential in the first unit cell to be much weaker 
than the bulk potential (c f Figs 3 and 4). The contri- 
butions to the a functions from excitations in different 
layers are shown in Fig 5 for photoemission from different 
layers. 

The mean free path is found to be considerably longer 
than obtained by Norman et al about 12 A rather 
than 3 A, at say 20 eV (Fig. 2). Measurements by the 
ITR-2PP technique ||| give a lifetime of r = 10 fa at 
an energy e = 3 eV above the Fermi surface. The mean 
free path is A = vt. Converting energy to velocity by 
mv 2 /2 = e gives a mean free path r = 103 A as com- 
pared to our result of about 17 A at that energy. This is 
an indication that our values rather are on the low side. 
It is however hard to know what is the correct conversion 
between energy and velocity at such low energies, which 
makes a comparison very uncertain. 

From Fig. 1 we see that the 2D losses occur only for 
small energies, at 5 eV the bulk losses take over. The 2D 
losses go to zero quite slowly, just like the bulk losses, 
but on another energy scale. If we only had 2D losses, 
the minimum mean free path would be long, about 20 A. 
The general behavior of the 3D mean free path follows 
a well known pattern. The mean free path has a mini- 
mum of about 5 A at an energy of 3-4 times the energy 
where the loss function has its center of gravity. We 
have used the Born approximation to evaluate the mean 
free paths. This may seem a very crude approximation 
at low energies. However the Born scattering expression 
with a basis of Bloch waves and Bloch energies rather 
than plane waves and free electron energies agrees with 
the GW approximation, which is commonly used also at 
low energies. Further it was shown by Campillo et al J^j] 
that plane waves and free electron energies was not that 
bad, as long as the energies in the dielectric function are 
well approximated. 

Our main concern is the behavior of the effective 
broadening function at small energies where it is dom- 
inated by the 2D losses. The 3D contributions set in at 
somewhat higher energies, and give a rather structure- 
less contribution. What we here for convenience call 2D 
losses is of course actually also a 3D effect since it comes 
form excitations of a coupled set of 2D layers. To allow 
a qualitative discussion we represent the ct2D functions 
by a rectangular distribution. Looking at Fig. 5 this 



may seem rather crude, but Fig. 8 shows that the cor- 
responding Pk (to) functions are not too different. The 
rectangular distribution allows an analytic solution (Eq. 
|l9| ) valid out to the cut-off uiq (ujq ~ 0.1 au ~ 3 eV). 
Pk (w) has only a fairly small tail beyond c^o ■ In Fig. 
9 we plot the total Pk (oj) function (sum over the four 
first layers, properly mean free path weighted), calcu- 
lated with the rectangular approximation and broadened 
with Lorentzians of different widths. We note the marked 
asymmetry. The asymmetry is described by an index 7 as , 
defined in Eq. ^0|. When (uj) derives from only one 
(rectangular) a-function "f as is a function of the singu- 
larity index olq, j as (ao). The index j as is then indepen- 
dent of the amount of Lorcntzian broadening T. If the 
Jib function has a power law singularity with singularity 
index ai, the asymmetry index contains the sum of the 
two indices, 7 as (ao + a£). 

In Fig. 7 we plot contributions to the loss function 
Pfc (w) from different layers. It is interesting that the 
first two layers give about the same contribution, while 
the contributions from the next two are tiny. In the left 
part of Fig. 6 we show the importance of the intrinsic 
contributions to olid- The behavior here is thus similar 
to what was found previously for metals. |^| In the right 
part of Fig. 6 we show the approach to the bulk limit 
curves. This approach is very slow for olid while, like in 
metals, it is fast for a^o- The slow approach for 012D of 
course comes from the slow approach to the bulk limit of 
the 2D fluctuation potentials (Fig. 4). 

In a paper by Liu, Anderson and Allen from 1991 [p7| , 
they discussed the lineshapes of Bi2Sr2BaCu20% along 
the r — X direction obtained by Olsen et al [^8| for 22 eV 
photons. They concluded that neither the Fermi liquid 
nor the marginal Fermi liquid theories could fit the slow 
fall-off of the spectrum at higher energies. Our results 
offer a possibility that the slow falloff may be due to in- 
trinsic creation of acoustic plasmons in a coupled set of 
CuO layers, an effect not present if only one CuO layer 
is considered. This broadening is mostly intrinsic, i e if 
we treat a 3D system we have an almost intrinsic effect. 
However most theoretical discussions concern an isolated 
2D system, compared to which we find an appreciable 
extra broadening from the coupling between the layers. 

The PES spectra change strongly when we go to the 
superconducting state. The main peak sharpens and a 
peak-dip- hump structure develops. This effect has been 
interpreted as a coupling of the 2D state to the (it, ir) 
collective mode. |lj| Here we find that this effect also 
can arise from the gapping of the loss function caused by 
the lack of low energy excitations in a superconductor as 
shown in Fig. 10. Without a more accurate model we 
find it difficult to decide which is the correct explanation, 
possibly it could be a combination of both mechanisms. 
Since the gapping of the loss function is related to the su- 
perconducting gap, also with our mechanism the hump 
will scale with the gap. It is clear that the experimental 
peak-dip-hump structure rides on a background which is 
not predicted by our expressions, nor by anyone else's. 
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Our theory is however rather schematic with its strict 
separation of a 2D and a 3D part, while in reality the 
bands are hybridized. If we extend our approach to a 
more detailed treatment of the underlying bandstructure, 
the background could well be strongly changed. Such an 
extension represents a very large numerical task but with 
the present pilot treatment we can at least start think- 
ing seriously about the difficult background problems in 
photoemission. 

In recent papers Joynt et al Q discussed a broaden- 
ing mechanism due to the interaction between the photo 
electron when outside the solid and the electrons in the 
solid. This is a different mechanism than in this paper, 
which adds additional broadening. Their discussion only 
involved the energy loss part and not the elastic contri- 
bution and can thus not be directly compared to exper- 
iment. We hence find their claims regarding pseudogaps 
uncertain. 

It should be stressed that we cannot claim any high 
quantitative accuracy. We have put in dispersion in the 
loss function using a crude approximation. Since how- 
ever dispersion is very important we think our predictions 
are substantially better than if dispersion had been ne- 
glected. We have only considered normal emission where 
the electrons come from the T point, while the interest- 
ing experiments concern electrons from the Fermi surface. 
However there is no reason that the effective loss func- 
tion should change qualitatively when we go away from 
normal emission. The behavior of the loss function when 
oj — > has been disputed. Most authors seem to believe 
the approach is linear, but there are also claims that it 
should be quadratic. |2^] If it were quadratic, the cor- 
responding a-function would start linearly rather than 
with a constant. However a (lo) would have to rise very 
fast to reproduce the behavior of the loss function for 
the (quite small) energies where it is known to be ap- 
proximately linear. Thus the pure power law behavior of 
Pk (uj) would be lost, but Lorentzian broadened curves 
would probably not differ much. Our fluctuation poten- 
tials are obtained by phase-shifting bulk potentials to 
make them zero at the surface, and define them as zero 
outside the solid. This procedure turned out to be fairly 
good in the metallic case, where we could check with 
more accurately calculated fluctuation potentials. Again 
this approximation is crude, but we believe it to be fun- 
damentally better than if we had used a step function on 
the bulk potential. Since the phase of the bulk potential 
is arbitrary, such a procedure would anyhow have been 
arbitrary. To calculate more accurate potentials is a very 
large numerical undertaking. 

One may also question the use of a bulk expression to 
estimate of the mean free path at the fairly low energies 
that we are concerned with, after all we found strong ef- 
fects when modifying the fluctuation potentials for sur- 
face effects. It does not seem easy to make a strong 
statement here, and we can only refer to "the state of 
the art" , that bulk mean free paths are successfully used 
in LEED and also in low energy life time calculations 



which are compared with time-resolved two-photon PES 
(TR-2PPE) experiments. |(| 
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V. APPENDIX A. DERIVATION OF THE PHOTO 
CURRENT EXPRESSION 

We will here derive Eq. g. The 2D and 3D parts in 
Eq. H factor, 

r (k,s h s 2 ) = J2 ( N 2D - 1, S2 \c t \N 2D ) t 3D (k,s u i) , 



where 



r 3D (k, Sl ,z) 
l + V 



E(N* B , Sl )+s k -H QB -h-V 



\N B )A\i). 



We note that (N 2B ~ 1, s 2 \ci\ N 2B ) is the basic part 
in the spectral function for the 2D system, and that 
the 2D and 3D parts are entangled through the index 
i. We have used Eq. ^| to eliminate the index s 2 , 
E - E (N 2D - 1, s 2 ) = E (N* B , Sl ) + e k . 
Expanding to first order in V we have |3lf] , 



V,f 



r 3D (k„s,z) =-I*(k|A|t 



(22) 



+ k 



V s 



e k -h 



where X is the self energy coming from a summation to 
infinite order in V, and we have used the relations 



(N B ,s\N B ) 



V, 



E(N B ,s) 



(N B ,s\V\N B } = V 8 (r). 

The energy argument of £ is uj s + e k . 

These results are only meaningful when u> s ^ 0. For 
the moment we take the excitation spectrum to have a 
gap, uj s > cjq, for all s except s = 0. For s = 0we have 



\(N* B AN B )\ 2 
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Also the s / terms in Eq. |22| have the exp (—a/2) 
factor, when we go beyond first order in V. 

We use plane waves for the parallel components of 
the fluctuation potential, V s (r) = e iQR V (q z , Q, z). Ne- 
glecting the reflected component, we similarly write the 
photo electron wave function as ipu (r) = e lKR- i/'j, £) (K,z). 
We further replace ImE by — iT, and absorb ReS in h. 
The photo electron energy is £k = (K 2 + fc 2 ) /2. We can 
now simplify the last term in Eq. ^ (c f Ref. ||) to 
become, 



pw , 



K-Q|(^f| V(q z ,Q,z) 



K 2 /2 -t z ~ {Very st (z) - Vo) 



All) 



where 



— =u s +e k Vo+iFi, (23) 



Ti = -ImE (fci,fc 2 /2) , kl/2 



£k - V a . 



t z = 



1 d 2 
2d^' 



pw - Q| is a plane wave 2D function, we have ne- 
glected the variation of the crystal potential in the lat- 
eral directions, the inner potential Vo is some average of 
V cr y st , and S° (fc, w) is the electron gas self-energy. For 
the ID Green's function we have approximately (see Ap- 
pendix C), 



1 



K 2 /2 -t z - {V cr yst («) - Vo) 



AiP<( Z< )tP>( Z> ). 



Here i(j< and ip> are damped Blochfunctions, i/>< decreas- 
ing towards the surface, and ^ decreasing towards the 
inner of the crystal (the crystal is on the positive half 
of the z axis), z > — max(z,z'), z < = min(z,z'), and 
the coefficient A is roug hly A = (in) . In our calcu- 
lations we will use the simplest possible approximation 
■0k {z) = cxp (inz), and ^< (z) — exp (— inz). 

The z part of the photo electron wave function is 



with 



f = f - Vo + *r 2 , r 2 = - ImS° (fc 2 , fc 2 2 /2) 



fc 2 /2 



£k - V . 



(24) 



Im E (e) we should choose the eigenvalue in the quasi 
particle equation (t + V crys t (r) — e^) tpu (r) = 0. For 
an electron gas this gives V crys t = ReS (fc, £&) + <j) DP 
since Vh — 0. The work function is by definition 
cf) = e F + Re S° (fcp, £fc F ) + DP . Since ReS°(fc,e fe ) 
varies fairly slowly with k out to about 2kp, we can take 
Vo = (f> — Ef = — \4>\ — £f- When we leave the elec- 
tron gas a reasonable definition for Vo is Vo — — \4>\ — W, 
where W > is the bandwidth. The maximum kinetic 
energy the photo electron can have outside the solid is, 
by energy conservation, u v hot — \4>\i corresponding to the 
energy u> p hot + W inside the solid. In our calculations 
we have taken \<f)\ = W = 0.15 au. Errors in this choice 
have a minor effect, and the relative error decreases with 
increasing photon energy. All calculations are made for 
a photon energy of 1 au. 

We consider only forward propagation for the Green's 
function from the excited layer to the surface (c f p2|), 
and have, 



T (k, S ,i) 



Vu 

-± k A i 



1 



+ -/ ^ (z')V(q z ,Q,z')^<(z')dz' 



J°° dz"^> [z") J dRe- J ( K - Q > R A (R,z") 
xi)i(H,z"), 



(k,s,i) 



UJ S 

i r° 



(fc„K|A|<> 



dz^ (z)y( 9z ,Q,z)V< (*) 



x (k,K — Q |A| i) . 

If we approximate the ip functions with plane waves the 

dipole matrix elements (k z ,Vi | A| and (k, K — Q | A| i) 

depend on the position of the excited layer through the 

factors exp [ik z zo^j and exp(iKZo). We note that the 

electron lifetime is r = 1/ (2r) and from Eq. ||we have 

2Imk z ~ 2r 2 / k z = 1/A where A = vr is the mean free 

path. Neglecting the recoil momentum Q picked up by 
the quasi-boson "s" the entanglement between the 2D 
and 3D parts disappears and we have the intuitively ex- 
pected result, 



Jk(z ,ujphot) = / J k (zo,uj)Pk(z 0l uj phot -uj)duj, 



(25) 



with 



We note that Ti and r 2 are different. 

We give a few comments on the relation between the 
electron energy inside the solid and outside. V crys t (z) 

is defined as V cryst = V H + Re £ (u) + <p DP , where Jl° {zq,^) =22 



b DP is the dipole contribution to the work function 
6, with <p defined as negative. For the argument e in 



J2( n id - l,s 2 \ci\N 2D ) (k,K\A\i) 



x5 {uj + E S2 - £k) 
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Vi 



i(^K — k z )z 



IK 



• I dze 
/o 



(k z — k) ; 



V(q z ,Q,z) 




innocent. Interlayer coupling in x° is absent in any one 
electron theory with a local potential, and thus e g in 
RPA. It is also absent in the static (ui = 0) case since 
this case can be described by DFT where the potential 
is local. 

We Fourier transform with respect to R, and separate 
into contributions from different layers 



where s — (q z ,Q). We have included the s — term 
and the common factor exp (— a). The photo current 
Jk(zo,U!phot) thus is a convolution between the sudden 
approximation 2D current J^ D (z ,ui), and an effective 
loss function Pk (zq, lu). 



x(z,z';Q) = 22xnn'[z,z';Q), 



(27) 



VI. APPENDIX B. DIELECTRIC RESPONSE 

Dielectric response is usually treated in the random 
phase approximation (RPA), and RPA has indeed proved 
extremely useful in many cases. | j33| For e g high T c ma- 
terials RPA may however not be good enough, and we 
will derive formal expressions without resorting to RPA. 
These expressions allow us to connect the energy loss re- 
sults to the screened potentials needed to discuss photo 
emission. The energy loss data are then taken from ex- 
periment. Some of our results can be found in Griffin's 
classic paper Q|, but not those which are crucial to our 
treatment. 

The response functions x , X, an d £_1 are defined from 
(in a schematic notation) 

pind ^Oytot ^yrext ytot yp^ nd -f- \f ex t ^—lyext 

This leads to the relations 

e -1 = 1 + vx, X = x" + X° V X- 

Since p md and V ext are exactly defined, no approxima- 
tions are involved in the definitions of x°, X> an d e _1 - 

We now specialize to two layers per unit cell. We 
choose the origin of the z coordinate at the center of 
the cell such that we have two layers at z = ±d. Wc 
write the response functions as, 

±i 

X°(r, r') = w(z— cm— dn)w(z'— cm— dn)x°(R— R') 



±i 

x(r, r') = w(z — cm — dn)w(z' — cm' — dn') 

mm' nn' 

><Xnn'(m - to'; R - R') 



We have assumed translational invariance in the layers, 
and that there are no transverse excitations, i e that 
the electrons always stay in the lowest transverse state 
(f>o (z) , w (z) = \4> (z)\ 2 . We have taken the overlap 
between w (z) and w (z + c) as zero, and neglected in- 
terlayer coupling in x° ■ This latter neglect is probably 



where 

Xn [z, z'; Q) = w(z — cm — dn)w(z' — cm — dn)x°(Q), 

m 

Xnn' (z, z'; Q) = w(z — cm — dn)w(z' — cm' — dn') 

mm' 

XXnn'(m - to'; Q). 

The integral equation x — X° + X° V X can be written as 
(suppressing the Q variable), 



Xnn'(z,z') = X n {z,z)5,, 



(28) 



/ Xn(z, Zl)v (Zl, Z 2 ) Xn"n'(z2, z')dz 1 dz 2 . 



This is the same result as in Eq. 1 in Griffin's paper |54| . 

We Fourier transform x°(z,z') with respect to z and 
z' , using discrete q z values and the orthonormal set 
{L- 1 /' 2 exp(iq z z)} 

X° n (q Z) q' z ) = \w (q z ) w (q' z )x°e<^) dn , 

where it is understood that q z and q' z differ by a recipro- 
cal lattice vector G, L is the length of the sample, and c 
the lattice constant. Similarly we have for Xnn' (qz, q' z ) 



Xnn'(q Z M = -W(q z )w(q' z ) Xnn' (& ) ^ ""e"^ 



(29) 



where 



We note that Xnn' (q z ) is a periodic function in g z , 
Xnn' (q z ) = Xnn' (q z + G) . Eqs. |7| and |9| give 

X (q z ,q' z ) = -w (q z ) w (q' z ) £ Xnn' (q z ) e^ dn e~< dn ' . 



(30) 
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We can separate out the w factors in Eq. g8| to obtain 
an equation for Xnn' (<?*) 



Xnn' (qz) = Xnn' + *£ m ^nm 2 (Oz) Xn 2 n> (<fe) 



(31) 



where Xnn' = X $nri and V is a 2x2 matrix periodic in 

Qz, 

Vnn' {q z ) = Y, v (& + G ) w ' (9- + G ) e- l ^ +G)d ("-"'). 



Eq. |3l] gives the matrix solution 

1 



Xnn' {qz) 



X 



l-V{q z )x°/c 



We write 



V ( n \-( V ^^)^ (qz) 
Vnn ' {qz) -\V 1 {qz)e-^) V {q z ) 



where 



V {q z ) = V h i {qz) = J% (?, + G) w2 (?* + G ) > 



and have 



~ r „ n _ 1 ( X1+X2 (xi -X 2 )exp(i<; 

Xnn>{q Z )- 2 I (xi _ X2)exp( _ #) Xl+X2 



Xi 



x° 



1 - (x°/c) (Vo + ^i) ' X2 1 - (X°/C) (Vb - Vi) ■ 



Here xi and X2 are functions of g z , Q and lo. From Eq. 
p0| the energy loss function becomes, 

v(q z ,Q)Imx{qz,qz;Q,u) 
_ ^fe,Q)^ 2 (gz) 

c 

x Im [Xl + X2 + (Xi - X2) cos (2g z d + {q z , Q))] , 

to be compared with the screened potential Im W = 
v (Imx) v. From Eqs. ^ and 29 we have 



ImW {Q,uj;z,z') 

±1 1 r 7r/c 



(32) 



7T- / V^(Q,g z ,z-dn) 

nn' / 



Xrm'(Q' W >9z 



xV* {Q,q z ,z' - dn)dq z , 



V{Q,q z ,z) 



G 



(33) 



There is thus no simple relation between IvaW {z, z') 
and the loss function unless the non-diagonal elements 
in Xnn'{qz) can be neglected. For typical values of q z 
and Q it however turns out that Vt/Vo is 0.2-0.3. Tak- 
ing xi = X2 and using the symmetries V{q z ,—z) = 
V{-q Zl z) = V*{q z ,z) and Xn{q z ) = Xn(-qz) we can 
write Eq. as 



2 ±1 



ImW {z, z') =J212~ I ' Vi{Qz,z-dn) 

i=l n ^° 

Imxofe) T/ , , , 

x Vi{q z ,z ~dn)dq z , 

c 

where V\ is the real and V% the imaginary part of V (q z , z) 
and xo = (xi +X2) /2. The real and imaginary parts 
turn out to give equal contributions to Im W. 

So far we have results for a set of coupled layers sitting 
in vacuum. We can take account of the embedding elec- 
trons (the 3D bulk excitations in our parametrization) 
by using 

X°{z,z':Q) =xg(z-z';Q) (34) 
+ w{z — cm)w{z' — cm)x2 (Q) • 

1 a 

This leads to a 3D bulk contribution in ImW {z, z'), and 
to screening of the 2D susceptibility x- The same screen- 
ing however appears also in the loss function, so we can 
forget about it in our problem. We note that the bulk 
screened potential can be anisotropic since Xb (<Zz> Q) can 



depend on both q z and Q, and not only on q 



Q 



We derived the relation between x and x° by solv- 
ing the integral equation x = X° + X° V X- This equa- 
tion can be written as x — X2D + X2dv?,dXi where 
X2D = X° + X° V 2DX2D- Here v 3 d contains no intralayer 
parts, while V2D only has intralayer contributions. Since 
X2D is available from many sophisticated theoretical cal- 
culations, it is interesting to have the relation between 
X° = X2D an( i X2D- We write 

X2D {z, z') =w{z)w {z') X2D, 



X2D { z , z ') =W{Z)W {z) X 2 D- 

The equation X2D = X° + X%.dX2D gives, 
X2D = x\d + X2D W ooX2D, where Woo = 
J w (zi) v {z\ — Z2] Q) w {Z2) dz\dz2- The desired relation 
is, 



X2D 



X2D 



1 + WooX2D 



where 
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VII. APPENDIX C. ON GREEN'S FUNCTIONS 
WHEN THE POTENTIAL IS COMPLEX 

Green's function theory is usually developed using 
a real potential. Here we will shortly summarize the 
changes from having a complex potential. With a con- 
stant complex potential the equation (Im/t > 0), 

\ + J?) G ^ Z ' Z ' ] ^ = S<yZ ~ 

has the solution (as is easily verified by direct substitu- 
tion), 



G(z, z'; k) = — e 
in 



IK\Z — Z 



With k a function of z the solution has the form, 

G (z, z'; k) = wg- (z<) g+ (z>) , 
where w is a constant (c f e g Arfken, |35|]), 

2 

w 



9- {z)g' + 0) - g'_ (z)g + (z) 



The boundary conditions are g- (z) — > for z 
and g+ (z) — > for z — > oo. 

In a slightly more general situation 



-co 



k{z) 



K\, z < 



Iiiik, > 0, 



K2, Z > 

we have 

<?_ (z) = e - lKlZ 6> (-z) + (a_e iK2Z + 6_e" lK2Z ) (z) , 
g+ (z) = e lK2Z (z) + (a+e iK12 + 6 +e - JKlZ ) 9 (-z) , 



«2 



For high energies ki ~ K2 and the results reduce to those 
of the first model with K=constant. 
Finally we consider a model with, 

, , f Kl, Z < T _ 

K ^ = \ K2 (z), z>0 ' Im ^> ' 



9- (z) - e-' iK1 ^ (-z) + (o_^> (z) + b^< (z)) (z) ; 



g + (z) = ^> (z) 9 (z) + (a+e"^ + b + e~^ z ) 9 (-z) . 

To see the difference between ip< (z) and -0k ( z ) we use 
the WKB approximation for the case z > 0, 

/> / \ i I K2(z')dz' /< / \ —i I K,2{z')dz' 



Since K2 (z) is complex, the two solutions are basically 
different, not just complex conjugates. For higher en- 
ergies and z > we have g_ (z) ~ ip< (z) , g + (z) = 
v> (z). 
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